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PREFACE 


It  is  my  policy  to  teach  nuclear  reactor  physics  as  a 
problem  oriented  course.  Such  problems  generally  require, 
as  input,  group  cross  sections  and  cross  section  related 
constants  such  as  the  diffusion  coefficient  or  neutron  age. 
However,  these  problems,  unlike  their  counterparts  in  the 
real  world,  usually  involve  only  one  or  a  few  energy  groups. 
Ibis  few  group  distinction  has  always  made  the  acquisition 
of  group  cross  sections  a  particularly  frustrating  experience. 
One  borrows  experimental  value  from  here,  a  calcu'J#£ed 
value  from  there,  and  combines  thenqwith  some  fast  cross 
sections  collapsed  by  hand  from  a  many  group  set  (■rfiich 
inevitably  doesn't  contain  all  the  elements  of  interest). 
Although  sophisticated  methods  are  always  available,  data 
for  one  problem  of  a  10  problem  assignment  rarely  warrants 
the  effort  to  obtain  themj  nor  does  the  accuracy  required 
in  a  school  problem  warrant  their  sophistication. 

The  code  described  in  this  report  is  an  attempt  to 
provide  reasonably  good  cross  sections  from  a  single,  easily 
used,  and  -unsophisticated  (therefore  inexpensive)  source.  Hie 
code  is  written  for  the  IBM  1620,  a  digital  computer,  installed 
in  11*6  (55/S)  of  the  268  schools  listed  in  ths  1961*  Oak  Ridge 
Institute  of  Nuclear  Studies  report,  SJJCAT10NAL  PRO  (RAMS  AND 
FACILITIES  IN  NUCLEAR  SCIENCE  AMD  ENGINEERING. 

ii 


Robert  H.  Hansen  carried  out  the  moments  reduction  study 
as  his  M.S.  thesis.  Ernest  Park  Sims  added  the  thermal  portion 
of  the  code  in  conjunction  with  his  M.S.  thesis  research,  and  put 
the  code  in  essentially  its  present  form.  In  addition,  acknow¬ 
ledgment  is  due  Siillip  A.  Terry  for  the  inclusion  of  fission 
cross  sections. 

Finally  I  feel  compelled  to  some  comment  about 
the  name  OLD  BARNYARD.  It  was  originally  coined  as  a 
joke;  however,  I  find  it  is  a  name  tdiich  the  student  remembers 
and  more  importantly  associates  with  cross  sections  and 
age.  I  have  yet  to  have  a  student  ask  "What  was  the  name  of 
that  code  which  generates  cross  sections?" 

Charles  J.  Bridgman 

Assoc  Prof  of  Nuclear  Engineering 

Air  Force  Institute  of  Technology 


Typed  by  Miss  Shari  L.  Johnston 


OLD  BARNYARD:  A  CROSS  SECTION  CODE  FOR  SCHOOL  USE 

ABSTRACT 

This  interim  report  describes  a  digital  computer  code 
■which  calculates  few  group  fast  and  thermal  nuclear  cross 
sections  and  constants.  It  performs  these  calculations 
using  the  moments  method  for  fast  flux  and  Wigner- Wilkins 
method  for  thermal  flux.  The  code  is  designed  for  use  on 
the  lift  1620  computer  and  is  intended  primarily  for  use  by 
students  and  professors  in  support  of  classroom  assignments. 

* 

INTRODUCTION 

Reactor  physics  instructors  who  assign  their  students 
criticality  or  associated  problems  are  continually  faced  with 
the  problem  of  providing  group  cross  sections  and  cross  section 
related  constants  such  as  diffusion  coefficient  and  age.  The 
alternative  to  providing  such  constants  xs  to  require  each 
student  to  find  or  calculate  his  own  values.  This  alternative 
is  a  guarantee  of  considerable  effort,  of  doubtful  educational 
value,  on  the  part  of  the  students  and  of  considerable  variation 
of  final  answers  for  the  instructor  to  evaluate.  Standard 
practice  in  most  first  year  courses  in  reactor  physics  is  to  use 
a  Maxwell  Boltzmann  distributicxi  to  determine  thermal  group 
constants  plus  experimental  values  as  available,  and  to  borrow 
epithermal  and  fast  group  values  from  some  published  set  such 
as  those  in  ANL  58001 .  Occasionally  an  instructor  will  have 
the  resources  to  obtain  fast  group  cross  sections  from  a 
moments  code  such  as  GAM2 .  However,  since  the  instructor 
usually  wants  to  use  only  a  few  groups  for  learning  purposes,  this 
last  option  is  inefficient  in  time  and  dollars.  More  likely  he  will 
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set. 

This  use  of  the  Maxwell  Boltzmann  distribution  for  thermal  cross 
sections  is  usually  adequate,  but  the  crude  collapsing  of  a  many 
group  set  for  fast  cross  sections  leaves  something  to  be  desired. 

As  fast  reactor  calculations  become  more  and  more  important  in  tne 
classroom,,  this  problem  is  amplified.  Further  when  cross  sections 
are  specified  in  this  manner,  one  can  almost  guarantee  a  mismatch 
in  values  such  as  La  /  D/E, .  All  of  this  uncertainty  is  very 
disconcerting  to  the  beginning  student. 

The  computer  code  described  in  this  report  is  the  r  isult  of  an 
attempt  to  provide  a  simple  easily  used  code  to  generate  few 
group  cross  sections  and  constants  for  school  use.  The  code 
calculates  both  fast  and  thermal  cross  sections  and  constants 
such  as  age  and  diffusion  length.  It  is  designed  for  execution  on  the 
IEH  1620  computer  idle  re  it  requires  about  lf>  minutes  for  execution. 

The  cross  sections  are  found  by  calculating  the  energy  dependent 
flux  in  the  infinite  homogeneous  mixture  of  the  isotopes  or  compounds 
specified.  The  input  cross  section  data  is  obtained  from  an  eleven 
group  library  of  fast  cross  sections  plus  2200  Vsec  thermal  values. 

The  energy  dependent  flux  is  calculated  by  the  neutron  transport,  moments 
method  for  above  thermal  energies  and  by  a  WLgne r-Wilkins  calculation 
for  thermal  energies  . 

The  moments  calculations  based  on  an  11  group  input  was  arrived 
at  after  a  systematic  reduction  in  the  fine,  group  structure*  until 
on  acceptable  compromise  between  accuracy  and  amount  of  calculation 


2. 


was  found  at  11  fine  groups.  The  above-thermal  calculations  were 
further  simplified  by  computing  resonance  absorption  with  a  resonance 
integral  taken  from  the  formula  given  by  Murray3  and  other  texts.  The 
Wigner-Wilkins  spectra  is  normalized  to  the  moments  spectra  at 
1.125>  ev.  The  details  of  these  calculations  are  given  in  section 
titled  "Theory." 

The  cods  itself  is  described  in  the  section  titled  "Hie  Code? 
along  with  simplified  operating  instructions  and  output  samples. 

Again  the  unique  feature  of  this  code  is  that  unlike  the  moment 
codes  and  thermal  codes  now  available  it  can  be  executed  on  a 
small  machine  (IEH  1620)  in  about  15  minutes  or  roughly  one 
dollar  of  computer  time. 


THEORY 


In  this  section  the  methods  of  solution  used  for  various 
important  calculations  in  the  code  are  presented.  Derivations  are 
cursory  and  less  important  or  routine  calculations  are  omitted.  An 
attempt  is  made,  however,  to  cite  appropriate  references  where  the 
interested  reader  may  find  a  detailed  discussion  of  the  subject. 
Epithermal  Age  and  Flux 

The  calculation  of  epithermal  age  and  f.ux  is  made  using  the 
moments  method;  i.e.,  the  moments  of  the  Fourier  transformed  Pa 
approximation  to  the  Boltzmann  neutron  transport  equation. 

The  derivation  of  the  equations  involved  is  straightforward 
but  long,  R.  H.  Hansen4  has  treated  them  in  great  detail.  The 
resulting  finite-difference  equations  to  be  programmed  are 
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and  N  =  total  nuclei  density,  nuclei  cm-3  x  10  34 

=  microscopic  removal  cross  section  for  the  n^  group,  bams 

aTn  =  microscopic  total  cross  section  for  the  n^1  group,  bams 

a°  .  n  "  zeroth  harmonic  of  the  microscopic  elastic  scattering 
'  cross  section  from  group  j  to  group  n,  bams 

a1  .  =  first  harmonic  of  the  microscopic  elastic  scattering  cross 

s  J  n  section  from  group  j  to  group  n,  bams 

r.  pth  derivative  of  the  m^  moment  of  the  Fourier  transformed 
mn  r 

flux  in  group  n 

•rn  =  neutron  age  to  the  lowest  energy  of  the  group  n,  cm3 

Xjj  =  fission  source  in  group  n 

In  addition  to  calculating  age, the  above  equations  generate  relative 
group  fluxes  through  the  These  fluxes  and  age  values  are  calculated 

down  to  lethargy  17  (O.UlU  ev).  All  values  of  0(u)  are  normalized  to 
0(16)  =  1  in  order  to  Join  with  the  thermal  flux  at  the  same  lethargy 
value. 

Resonance  Escape  Probability 

The  calculation  of  resonance  escape  probability  follows  the  treatment 
given  in  Murray  1  s  'text5 .  The  scattering  per  resonance  atom  is  calculated 
and  the  value  obtained  is  used  to  select  an  empirical  expression  for  an 
effective  resonance  integral.  The  empirical  equations  used  are  taken  from 
Isbin6  and  are  reproduced  below: 
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where  £g  =  macroscopic  scattering  cross  section  for  the  mixture  over 
the  resonance  region,  cm-1 

-  nuclei  density  of  nuclide  A,  cm-3 

The  value  of  the  resonance  escape  probability  is  then  found  from 

P  -  exp  r-(RI)eff  NV5  2  s3 

where  §  Es  is  for  the  mixture. 

This  value  of  p  is  then  introduced  into  the  equations  (l)  by  multi¬ 
plying  a^a  by  p  (assonance  takes  place  in  group  8).  This  a.-tifioe  has 

the  effect  of  reducing  the  downs cat ter  to  group  9  by  appropriate  factor 
without  violating  neutron  conservation  within  the  equations. 

Should  the  user  so  desire,  he  may  input  his  onfc.  value  of  (Kl)e^ 

to  be  used  in  the  code. 

Absorption  Parameter 

The  absorption  parameter  is  defined  for  a  single  nuclide  by 
2A7a(B)/  J/kT 

4  '  — -  ft) 

vdiere  A  *  mass  of  the  nuclide,  amu 

oa(E)  =  microscopic  absorption  cross  section  at  energy  E,  bams 


cr  (E)  *  microscopic  scattering  cross  section  at  energy  E,  bams 
s 

k  =  Boltzmann’s  constant 
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The  equation  i c  solved  by  setting 

3 

N(x)  »  x3e-X  M(x) 


which  converts  it  to 

xM  1 '  (x)  -  (3-2x3  )  M*  (x)  -  AM(x)  -  0 
This  equation  is  than  solved  by  the  power  series  method. 
The  result  is 

CJ 
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In  the  code  the  value  of  N(x)  and  x3  N(x),  which  is  0(u),  are 
calculated  at  every  0.2$  lethargy  interval  beginning  at  u  =  23  and  con¬ 
tinuing  to  u  =  16.  The  calculation  employs  equation  (10)  and  adds  terms 
until  contributions  of  additional  terms  are  less  that  10 6  relative  value. 
All  thermal  fluxes  are  normalized  to  0(16)  ®  1  in  ordsr  to  match  the 
epithermal  fluxes  from  the  moments  calculation. 

Thermal  Constants 

The  thermal  range  in  the  code  is  arbitrarily  defined  as  the  interval 
0  to  1.125  ev  which  corresponds  to  <*  to  16  on  the  lethargy  scale.  ADI 
thermal  constants  are  calculated  for  this  interval. 
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Host  Probable  Neutron  Velocity.  The  most  probable  thermal  neutron 
velocity  is  found  by  locating  the  peak  of  the  N(x)  vs  x  curve  generated 
in  the  thermal  flux  calculation.  Since  values  of  N(x)  are  calculated 
at  intervals  of  0.25  lethargy  rather  than  continuously  the  value  of 
vp  found  is  not  exact. 

Average  Neutron  Velocity.  The  average  neutron  velocity  is  found 


from  the  formula 


J  x  N(x)dx 


_  o 

x  =  “eTTT 


J*  N(x)dx 


where  x  is  again  the  normalized  velocity  and  6.6?  is  that  value  of  x 
corresponding  to  lethargy  16.  The  integral  overN(x)  is  calculated 
analytically  by  integrating  the  Wilkins  equation.  The  integral  over 
xN(x)  is  found  by  the  trapezoidal  rule  during  the  thermal  flux  calcula¬ 
tion.  Once  x  is  found  the  most  probable  neutron  velocity  is  given  by 

Vp  *  x^l.65  x  10 4  T  ft 

where  T  is  again  system  temperature. 

Thermal  Diffusion  Length  Squared.  The  value  of  thermal  diffusion 
length  squared  is  calculated  from  the  equation  9 

a  Ss  * 

L2  -  - 1 -  (] 

3  E*  E 
T  ao 


°a  x 

L  3aao  C(Gs  +aao^>N]3 

where  L  =  thermal  diffusion  length,  cm 

and  where  a40  =  [xa»  (x)3x=x^  ,  Xp  is  most  probable  value. 
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Here  the  cross  sections  and  nuclei  density  are  input  data  to  the 


program  and  x  is  calculated  as  described  earlier. 

Migration  Area.  The  total  migration  area  for  a  neutron  is  found  by- 
adding  the  moments  method  age  to  lethargy  16  to  the  value  of  L?  for  the 
lethargy  interval  ®  to  16. 

Thermal  Cross  Sections.  The  2200  m/sec  cross  sections  for  appro¬ 
priate  nuclide 3  are  input  data  to  the  code.  These  values  are  adjusted 
by  the  code  to  account  for  deviation  of  the  thermal  spectrum  from  Max¬ 
wellian.  In  particular  the  proper  absorption  cross  section  for  the 
thermal  group  is 


and  the  transport  cross  section  is 

°tr  =  +  astou)  to) 

where  the  term  (1-13)  is  also  input  data. 

Cross  Section  Collapsing 

Cross  section  libraries  for  the  code  have  a  twelve  group  structure; 
eleven  epithermal  and  one  thermal.  Only  the  epithermal  groups  may  be 
collapsed.  The  thermal  group  is  always  left  intact  and  is  output  in 
the  same  manner  regardless  of  the  number  of  epithermal  groups. 

The  eleven  epithermal  fine  groups  may  be  collapsed  into  any  number 
of  broad  groups  from  one  to  eleven.  The  boundaries  of  the  broad  group 
are  restricted  to  those  values  which  are  boundaries  of  the  fine  groups. 

Broad  group  cross  sections  are  generated  from  the  formula 
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where  o =  microscopic  cross  section  of  the  n  kind  for  a  broad  group 
01  =  total  flux  in  the  ith  fine  group 

cr^  -  microscopic  cross  section  of  the  kind  for  the  ith  fine 
group 

An  exception  to  the  method  of  equation  (18)  is  the  calculation  of 
the  transport  cross  section.  Here  the  use  of  o^.r  in  a  reciprocal 

manner  through  the  diffusion  coefficient  is  anticipated  so  t^at 


is  used. 

A  number  of  different  cross  sections  are  output  by  the  code  and 
the  interrelationships  which  exist  among  them  is  not  always  obvious. 

The  following  explanation  of  more  obscure  relations  may  prove  useful 
to  the  user. 

Total  Transfer.  The  group  to  group  total  transfer  cross  section  is 

that  value  which  when  multiplied  by  the  losing  group  flux  gives  the  rate 

of  total  transfer  from  the  losing  group  to  the  gaining  group.  Thus  total 

transfer  from  j  to  k  consists  of  the  sum  of  inelastic,  elastic  (P  ),  and 

o 

twice  n-2n  from  j  to  k.  Note  that  elastic  (Pi)  which  is  jl  times  elastic 
(P0)  does  not  contribute  to  total  transfer. 


Removal.  The  removal  cross  section  is  applied  only  to  a  broad  group 
and  is  that  value  which  when  multiplied  by  its  broad  group  flux  gives  the 
rate  of  removal  of  neutrons  from  the  group.  Thus  the  removal  cross  sec¬ 
tion  is  the  sum  of  absorption,  inelastic  removal,  elastic  (P  )  removal, 

and  one  times  the  n-2n  removal  cross  sections.  Note  that  in-group  cross 
sections  are  not  included  as  they  do  not  remove  neutrons  from  the  group. 


■  i 
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THE  CODE 

Ihe  digital  computer  code  which  carries  out  the  calculations 
described  in  the  previous  section  is  written  in  the  Kingston  version 
of  Fortran  H10  (sometimes  called  "Kingstran").  The  code  consists  of 
two  chains  which  must  be  run  sequentially.  Ihe  first  chain  computes 
the  energy  dependent  flux,  both  fast  and  thermal.  The  output  of  this 
chain  consists  or  the  group  flux  for  each  of  the  11  fine  fast  groups 
plu$  9  thermal  values  at  half  lethargy  increments  as  well  as  age, 
diffusion  length  squared,  resonance  integral  parameters,  and  average 
and  most  probable  neutron  velocities  in  the  thermal  region.  In  add¬ 
ition  the  chain  one  output  includes  a  plot  of  the  flux  as  a  function 
of  lethargy  over  the  entire  range  or  interest,  ihe  second  chain 
collapses  the  cross  sections  to  the  desired  group  structure  using 
the  flux  calculated  in  chain  one.  ihe  operator  may  choose  between 
microscopic  and  macroscopic  cross  sections  for  chain  two  output. 

A  listing  of  the  source  decks  for  chains  one  and  two  is  contained 
in  Appendix  a.  The  source  decks  shown  there  have  been  executed  on 
an  IBM  1620  with  UOK  memory,  care  input-output,  floating  divide, 
indirect  address  and  the  usual  additional  instructions (TNS,TNF,  MF). 
Operating  Instructions 

1.  Set  all  console  switches  to  PROGRAM  except  the  Parity  switch 
is  set  to  STOP. 

2.  Place  special  Old  Barnyard  short  subroutines  in  tne  read 
hopper.  Press  RESET,  LOAD. 
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3.  When  last  card  of  subroutines  is  reached,  the  console  will  give 
READER  NO  FEED  lig^t.  Press  READER  STaRT  and  read  in  last  card. 

U.  Remove  subroutines  from  the  out  hopper.  Place  Old  Barnyard 
chain  one  object  deck  in  the  read  hopper.  Start  the  PUNCH.  Place 
data  cards  for  chain  one  on  top  of  chain  one  object  deck  in  read 
hopper.  Press  LOAD. 

f>.  Chain  one  program  will  run  automatically.  When  last  data  card 
eaters  reader,  console  will  give  READER  NO  FEED  ligfct.  Press  - 
READER  START  to  read  in  last  card. 

6.  Console  typewriter  will  indicate  the  beginning  of  each  section 
of  the  code.  At  the  end  of  the  thermal  calculation  the  typewriter 
will  type  "SET  SWITCH  1  ON  TO  PRINT  OUT  SPECTRUM.  PRESS  START. "  If 
a  plot  of  the  total  neutron  spectrum  is  desired  set  switch  1  ON 
and  press  START.  If  no  spectrum  plot  is  desired,  set  switch  1 
OFF  and  press  START. 

7.  At  the  end  of  chain  one,  the  console  typewriter  will  type 
"END  OF  CHAIN  1.  TO  COLLAPSE  CROSS  SECTIONS  LOAD  CHAIN  2." 

If  no  cross  sections  are  desired  program  is  finished.  If  cross 
sections  are  desired  place  Old  Barnyard  Chain  2  object  deck  in 
the  read  hopper.  Place  data  cards  for  Chain  two  on  top  of  it. 

Press  RESET,  K)AD.  Chain  two  will  run  automatically.  When 
last  data  card  enters  reader*  console  will  give  READER  NO  FEED 
ligjrt.  Press  READER  START  to  read  in  last  data  card.  Remove  Chain 
two  and  data  dards  from  out  hopper. 

8.  When  chain  two  is  finished  the  typewriter  will  write  "END  OF 
PROGRAM. "  Remove  output  cards  from  output  hopper.  It  is  not 
necessary  to  press  NON  PROCESS  RUN  OUT. 

1U 


9.  List  output  answer  cards  on  the  hO?  lister.  Put  switch 
number  k  at  the  right  end  of  the  machine  UP.  This  causes 
automatic  advance  of  a  fresh  sheet  of  paper  for  each  new 
reaction  of  the  output  answers.  If  the  paper  advances  at 
•unwanted  times  (as  it  sometimes  does)  put  switch  U  DOWN.  (ihe 
best  plan  is  to  get  switch  k  UP  for  chain  1  output,  DOWN  for 
chain  2  output). 

10.  To  run  a  second  problem  go  back  to  step  (l)  and  repeat. 
NOTE:  If  anything  should  go  wrong  during  a  run  (such  as  an 
erroneous  data  card,  or  a  I/O  error,  or  a  SKIP  CHECK)  the 
particular  chain  in  use  may  be  restarted  by  the  following 
prooedure: 

1)  Press  STOP 

2)  Remove  remaining  data  cards,  if  any,  from  the  read 
hopper. 

3)  Clear  the  reader  of  any  internal  data  cards  by 
pressing  NON-PROCESS  RUN  OUT. 

h)  Replace  corrected  data  cards  in  the  read  hopper. 

5)  Press  RESET,  INSERT. 

6)  T^ype  on  the  console  typewriter  the  numbers  h900936. 

7)  Press  RELEASE,  START  or  hit  the  R-S  key  on  the  type¬ 
writer. 

8)  Press  READER  START  to  read  in  data  cards 

9)  Program  will  run. 


The  following  data  cards  are  required  to  operate  OBY  Chain  one. 
Unnecessary  decimal  points  may  be  omitted  on  all  cards. 

Card  No.  1: 

This  card  contains  65  columns  of  alphameric  data  of  the  user’s 
choosing.  Identification,  problem  title,  date,  user's  name,  etc. 
are  possible  information  to  be  placed  on  this  card.  The  data  on 
this  card  will  be  typed  out  on  the  console  typewrxuer  and  will  be 
punched  at  the  head  of  the  output  data. 

EXAMPLE i  Checkout  of  OBI  VIII  for  Carbon.  25  May  66  Wetzler. 

Card  No.  2: 

This  card  carries  five  (5)  numbers  as  follows,  in  any  desired 
formats 

a)  A  number  between  two  (2)  and  twenty-one  (21)  indicating 
the  source  of  neutrons  which  the  user  wants  to  use.  The  meaning 
of  these  numbers  is  shown  below* 


Number 

2 

3 

h 

5 

6 

7 

8 
9 

10 

11-21 


Source 

235 

U_  fission  (Cranberg  Spectrum) 

Tr^;: -fission 

PlCf' 

P’S 

Cf2*2 

Pu-Be  (Whitraore-Baker) 

Pu-Be  (Cochrane-Henry) 

Ra-Be  (Rill) 

The  user’s  own  source  (See  Cara  No.  U) 

A  unit  source  in  group  1-11  respectively 

EXAMPLE:  15  =  a  unit  source  in  group  $. 
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b)  A  second  number  indicating  the  number  of  nuclides  to  be 

used  in  the  problem. 

EXAMPLE:  2  for  water  (H  and  0) 

2  for  BeO  (Be,0) 

2  for  G, Hg  (C,H) 

3  for  an  AlgOg-I^O  mixture  (A1,0,H) 

1  for  pure  carbon  (C) 

c)  A  number  giving  the  total  nuclei  density  for  the  problem 

in  atoms  per  cm3  times  10-34. 

EXAMPLE:  for  carbon  0.08023 

for  HgO  0.1003 

(Note:  See  ANL  380(5  page  30) 

d)  A  number  giving  the  temperature  of  the  system  in  'degrees 
Kelvin. 

e)  A  number  indicating  whether  or  not  a  resonance  escape 
calculation  should  be  performed;  as  below: 

1  «  Yes 

2  =  No 

Card  No.  3* 

This  "card"  is  a  small  deck  of  cards  containing  the  source 
data  Cor  the  program.  It  is  included  as  a  part  of  the  input 
library  and  normally  requires  no  preparation  on  the  part  of  the 
user.  A  listing  of  the  sources  in  the  authors'  library  is  shown  in 
Figure  1. 

Card  No. 

This  card  is  required  only  if  the  first  number  on  card  no.  2 
was  "10";  i.e.,  only  if  the  user  wants  to  read  in  and  use  his  own 
source.  The  card  will  contain  eleven  numbers  in  any  format 
representing  the  source  in  each  of  the  first  eleven  groups  of  the 


.5  1.  1,5  2.  2.5  3.  5.  8.  12.  16.  17. 

THE  SOURCE  USED  IS  U235  FISSION  <  CRANBER6  SPECTRUM). 

2  .023023  .10824  .21044  .23139  .18048  .11483  .12439  .0072094  .0  .0  .0 
THE  SOURCE  USED  IS  U233  FISSION. 

3  .021495  .10472  .20811  .23202  .18252  .11674  .1267  .0073845  .0  .0  .0 

THE  SOURCE  USED  IS  PU239  FISSION. 

4  .024567  .11047  .21062  .22996  .17917  .11412  .1239  .0071986  .0  .0  .0 
THE  SOURCE  USED  IS  PU241  FISSION. 

5  .026852  .11615  ,215  .22949  .11077  .1190  .0068574  .0  .0  .0 
THE  SOURCE  USED  IS  CF252  FISSION. 

6  .046142  .14786  .2292  .21936  .15472  ,09531  .099137  .0055827  .0  .0  .0 
THE  SOURCE  USED  IS  PO-3E  ( V.'HI  TMORE-BAKER )  . 

7  ,2342  .357  .261  .1025  .036  .0085  .0  .0  .0  .0  .0 
THE  SOURCE  USED  IS  PO-BE  ( COCHRAN E-HENRY ) • 

8  .257  t 369  .2514  ,079  .028  .010  .0  .0  .0  .0  .0 
THE  SOURCE  USED  IS  RA-BE  (HILL). 

9  .239  .365  .189  .109  .065  .033  .0  .0  .0  .0  .0 
THE  SOURCE  USED  IS  YOUR  OWN  AS  SHOWN  BELOW. 

10.  0.  G.  0.  0*  0.  0*  0.  0.  0.  0.  0. 

THE  SOI  RCE  USED.  IS  A  UNIT  SOURCE  IN  GROUP  NO. 

1C  .0  .0  .0  .0  .0  .0  .0  .0  .0  .0  .0 

Code* 

First  line:  Upper  lethargy  limits  of  the  11  fine  groups. 

Second  line:  Identification 

Third  line:  Identification  number  followed  by  the  percent  of  the  source 
found  in  the  11  fine  groups,  in  the  order  S!  j 

Figure  1:  Input  Library,  Sources 
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program  group  structure. 

Card  No.  5: 

This  card  is  necessary  only  if  the  fifth  number  on  card  no.  2 
was  lj  i.e.,  if  e  resonance  escape  calculation  is  to  be  performed. 
The  card  will  contain  three  (3)  numbers  as  follows:  (any  format) 

a)  A  number  indicating  whether  the  resonance  nuclide  present 

is  Th^^,  or  w  (Tungsten)  as  below: 

1  =  D  or  W 

2  =  Th 

1  =  any  mixture  of  the  above  3 

b)  A  number  giving  the  nuclei  density  x  10^*J  of  the 

cm 

resonance  nuclide(s)  in  the  problem. 

c)  A  number  indicating  whether  the  effective  resonance 
integral  is  being  supplied  by  the  user,  or  whether  it  is  to  be 
calculated  by  tile  program,  as  below: 

0  =  to  be  calculate 
non-zero  =  value  of  to  be  used. 

Card  No.  6s 

This  card  contains  five  pieces  of  data  as  follows: 

a)  1 $  columns  of  alphameric  identification  to  identify  the 
n*^  nuclide  of  the  problem.  This  is  optional,  for  the  user’s 
use  only,  and  may  be  left  blank. 

xt. 

b)  &  number  giving  the  nuclei  density  of  the  n™  nuclide 

in  1  ftfarcS  x  1CT34J  in  the  problem. 
cm° 


c)  A  number  indicating  whether  the  nuclei  is  a  moderator, 
fuel,  or  other,  as  below: 

1  =  moderator 

2  =  fuel 

3  M  other 

d)  A  number  indicating  whether  micro  -  or  macroscopic 
cross  sections  are  to  be  output  for  this  nuclide.  This 
number  is  necessary  only  if  cross  section  averaging  is  to  be 
performed.  (Chain  two  of  the  program).  The  number  should  be: 

1  =  micro 

2  *  macro 

^l)  A  number  indicating  whether  fission  cross  sections  are 
to  be  included: 

1  *=  yes  (fission) 

2  =  no  (fission) 

Card  No.  7» 

This  "card"  is  a  deck  of  cards  bearing  the  11  group  plus 
thermal  cross  sections  of  the  n^  nuclide.  This  deck  is  selected 
from  the  program  library  and  requires  no  preparation  on  the  part 
of  the  user.  Some  examples  of  the  author’s  library  are  shown  in 
Figures  2  and  3. 

Cards  6  and  7  are  required  for  each  nuclide  ip;  the  program. 
The  ordering  of  nuclides  is  completely  arbitrary.  Figure  k 
shows  the  arrangement  of  the  entire  input  deck  for  Chain  one. 

Notes  on  Chain  One  data: 

a)  The  total  nuclei  density  on  the  first  data  card  should 

be  the  sum  of  the  individual  densities  on  the  separate  card  no.6's. 

b)  The  nuclei  density  of  the  resonance  nuclide  on  card 
number  5  should  agree  with  that  on  card  no.  6  for  the  absorber 
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CARBON®  62®0® 

2.0000®'  6.0000  ® 

4.6774E-02  2.7711E-0 

8.0933E-03  4.6132E-0 

5748E-G1  1.1047 

5.868  1E-0 1  -5.6251E-0 
1.2910  1.0688 

8.17092-01  -5.6305E-0 
3.7440  1.1643 

2.4442E-01  -2.1432E-0 
A. 5242  9.4954E-0 

7 .4309E-01  -6.5822E-0 

Code: 

a)  Identification  of  isotope  (element).  1st  l£  columns. 

b)  Number  of  fast  cross  section  enteries  below  (beginning  on  2nd  line) 

c)  E^ithermnl-fast  absorption,  o  =  negligible;  1  =  not  negligible 

d)  Type  °f  inelastic  scatter:  1  =  inelastic ;  2  =  (n,2n);  3  =  both;  Ij.  «  neither. 

e)  Iherraal  average  logarithmic  energy  decrement,  c 

f)  Thermal,  total,  microscopic  cross  section^  0. 

g)  Thermal,  absorption,  microscopic  cross  section, 

h)  Thermal,  absorption  parameter,  A 

i)  Thermal,  scattering  microscopic  cross  section,  0 

j)  Thermal  (l-u)  where  a  =  2/3A 

k)  Thermal  value  of  the  product  of  neutrons  per  fission  and  the  microscopic, 
fission  cross  section,  v0f 

l)  Number  of  groups  from  which  neutrons  inelastically  scatter 

m)  Number  of  groups  to  which  neutrons  inelastically  scatter 

n)  Microscopic,  fast,  inelastic,  cross  sections  in  the  following  order: 


an 

a  5 

oin  ; 

_in 

...  0 

i^i 

i-t> 

1-.7 

in 

0  } 

ain  } 

„in 

•  •  •  0 

p 

1®  .158®  4, 
[o.oooo 

2  3.1662E- 

3  5.5135E- 
3.2908E- 

1  1.C714 

6.5028E- 
1  2.1883 

3.2552 E- 
1  4.5213 

1  _1  •  8  5  7  7  E  - 


l]® 


.81®  .004®.  02® 

3.5757E-02 
■02  0.0000  ^ 

•03  O.OOOOj^ 

•01  -3.3102E-01 
9.6768E-01 
■01  -4.7099E-01 
1.8060 

•01  -2.4896E-01 
9.4884E-01 
•01  -1.6455E-01 


4.806®  .9444^  ®l 

1.0514E-01  7.8474E-02 

0.0000  4.0467E-03 

l.lOOOE+Ol'*  1.0000® 
9.8070E-01  1.3262 

7.1394E-01  --8.5887E-01 
1.7890  1.5656 

9.6579E-01  -7.0608E-01 
4.3996  9.8826E-01 

1.8565E-01  -1.6434E-01 
3.9669  1.4432 


o) 

P) 

q) 


w 

Microscopic,  P0  and  Pa , 


Number  of  groups  from  which  neutrons  elastically  scatter 
Number  of  groups  to  which  neutrons  elastically  scatter 

elastic  scatter,  cross  sections  in  the  following 

Po 

Oi-a 
Po 

°2-3 


order: 


Pr?i  J 


_  Pi 

ai-=i 


Po 

°a-a  > 


Pi 

<h-a 

Pi 

*7s-4 


Figure  2:  Input  Library,  Carbon  Cross  Sections 
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®  CS >  ©  <3?  ®  m 

•JRANIUM235  116  1  1  .008^704^ 

£l591E+0i  .1204E+01  .1303E+01 

.1598E+C1  .4071E+01  .2310E+02 

{• 8  E+01JS>.[6272E~03  •  2387E-01 

.4077  .4232  .1293E-01 

.3119  .6108  .4773 

.0  .2376E-C1  .1745 

. 3038E-01  .8512E-04  .0 
.4155  .8215  .4055E-01 

.0  .1156  .3104 

.1955E-07  .0  .0 

.2752E-03  .1633E-07  .0 
.7O98E-03  .2034E-C6  .0 
£3! 30E+01  . 1062E+02  .9810E-02 
• 2267E-01  -.4294E+01  .4619E+01 
.3860E401  .6376E+01  .3752E-01  - 
.4281E-01  -« 1068E+01  .4374E+01 
• 7348E+01  .3865E+01  .2801E-01  - 
•  2631E-01  -• 23  56E-0 1  .1048L+02 
. 1 241E+02  .1308  .2710E-C1 

.1  167  -.1038jv# 


689®  21582 
1305E+0 1 
7208E  t-02 
1738 
3459E-04 
5816 
4161 
0 

11  75E--03 

8789 

1965 

of  t 

• 1247E+01 

•1191E+02 

.9214 

.4665E+01 

.8345E-01 

.1077 

-.2493E-01 


P  15.^.99  72^1442.5® 

.  1255E+01  .1213E+01  m 
,8172E+02j^[7  E+0l]^ 

t4312  .5266 

.5427E-02  .8050E-01 
>  20 1 3E-G1  .  5497E-04 
.4959  .7476 

•6236E-01  .2587 
.6023E-08  .0 
.5527E-01  .1671E-03 
> 1032E  +  01  «  8628E-0 1 
6694  t..  1779 

11  E+02r^l  E+Ol] 


i© 


o  5203E+01 
.3605E-01 
. 3881E+01 
.5143E-01 
S 1043E+02 
•  2 195E-01 
. 1385E+02 


. I733E+Q2 
-.2 924E+01 
.5201E+01 
-.5365 
. 1182E+01 
-.1834E-01 
.2227 


[.5753E+01  • 3828E+01  .3786E+01  .3528E+01  .3247E+01  .3028E+01 
• 3728E+01  • 7827E+01  .3763E+02  .1012E+03  .17C1E+03]© 


Cede  i 
a)  -  q) 


r) 

s) 


t) 

u) 

v) 


Same  as  in  Figure  2 


t?. 


Microscopic,  fast,  absorption  cross  sections,  oi. 

Microscopic,  fast  vur  in  the  order  vo, ,  vu? 

Note*  These  cross  sections  must  not  be  included  for  chain  one 
but  must  be  included  for  chain  two. 

(Not  shown)  Number  of  groups  in  which  (n,2n)  occurs 
(Not  shown)  Number  of  groups  to  which  (n,2n)  scatter  occurs. 
(n,2n)  microscopic  cross  sections  in  same  order  as  inelastic 
Note:  t),  u),  v)  would  appear  in  that  order  just  before 
entry  q). 


Figure  3*  Input  Library,  U23f>  Cross  Sections 
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Same  for  nuclide  #3 


L  kU  -^^^Kuclei  densit 


for  nuclide  #2 


density  and  cross  sections, 
„  nuclide  #1 
►Resonance  data  (optional) 


►User’s  ovn  source  (optional) 

♦Source  deck  (supplied) 

“Source  number,  no.  of  nuclides,  total  nuclei  density,  T,  etc. 
♦Problem  ID,  title,  users  name,  date,  etc. 


FIGURE  4:  CHAIN  ONE  INPUT. 


or  -with  the  sum  of  those  on  several  card  no.  6’s  if  there  is  more 
than  one  absorbing  nuclide. 

Preparation  of  Input  Data  Ghain^Two 

Chain  Two  of  Old  Barnyard  requires  only  one  additional  data 
card.  This  card  bears  up  to  twelve  numbers.  The  first  is  the 
number  of  broad  groups  for  which  cross  sections  are  to  be 
calculated.  The  remaining  numbers  are  the  numbers  of  the 
lowest  fine  groups  in  the  selected  broad  groups.  Hie  last 
of  these  should  normally  be  eleven  indicating  that  the  last 
broad  group  extends  down  to  . UlU  ev.  Thus  if  the  first  number 
is  five  there  will  be  five  more  numbers  on  the  card.  If  the 
first  number  is  11,  there  will  be  11  more  numbers  on  the  card. 
Examples: 

5  2  U  6  8  11 

11123U56789  10  11 


This  additional  card  is  bhe  first  input  card  for  dhain 
two.  It  is  followed  by  the  collection  of  cards  6  and  7  from 
chain  one  as  shown  in  Figure  5>. 

Output 

A  sample  output  of  the  code  is  shown  in  Figure  6.  This 
output  is  for  a  homogenized  mixture  representing  an  MTR  type 
core.  Figure  6a  shows  the  output  of  chain  one.  Figure  6b 
shows  the  flux  plot.  Figures  6c,  d  and  e  are  the  chain  two 
output.  In  this  example,  the  user  has  asked  for  three  group 
macroscopic  cross  sections  which  have  boundries  at  fine  groups 


The  code  in  the  form  presented  in  Apendix  A  has  been  used 
routinely  by  several  first  year  graduate  students  during  the 
past  academic  term.  These  students  had  access  to  a  set  of 
instructions  essentially  equivalent  to  those  contained  here, • 
and  have  used  the  code  successfully  without  any  other  tutoring. 

The  code  executes  in  about  1$  minutes  on  the  authors’  IM  1620. 
At  least  five  minutes  of  this  time  is  involved  in  computing  the 
flux  plot.  However,  this  plot  feature  is  believed  to  be  very 
worthwhile  from  an  educational  point  of  view.  Actually  seeing 
the  energy  dependent  flux  for  a  particular  problem  affords 
the  student  considerable  insight  into  the  physics  of  that 
problem. 

Status 

There  are  several  minor  bugs  in  the  program  as  presented  in 
Appendix  A.  The  most  annoying  of  these  is  failure  of  the  plot- 
routine  for  a  fast  mixture  (say  pure  uranium).  Another  dis¬ 
concerting  feature  is  that 

otf  -  a,+  aia  +  aP0  -  crPx 

This  is  probably  due  to  the  reciprocal  formulation  of  atr  • 

These  problems  are  being  corrected  now  and  in  addition  better 
temperature  and  most  probable  thermal  velocity  calculations  are 
being  incorporated.  Also,  it  is  planned  to  incoiporate  simple 
fast  fission  factor  and  helgogeneous  cell  calculations. 


Ihe  final  code  along  with  a  comprehensive  theory  and 
operating  instructions  will  be  published  as  the  next  and  final 
report. 


NETF  CORE  COMPOSITION  FOR  DB 


THE  SOURCE  USED  IS  U235  FISSION  (CRANBERG  SPECTRUM ) « 
SYSTEM  TEMPERATURE  IS  293.0  KELVIN 
ATOM  FRACTIONS  ARE 


URANIUM235 

.00137621 

URANIUM  38 

.00015252 

HYDROGEN 

.48262000 

OXYGEN 

.24131000 

ALUMINUM 

.27583353 

GROUP 

LETHARGY 

ENERGY*  EV 

SOURCE 

RELATIVE  GROUP 
FLUX*  /UNIT  U 

1 

.50 

6.0653E+06 

2.3023E-02 

5.4905E-01 

2 

1.00 

3.6788E+06 

1.0824E-01 

2.2786E+00 

3 

1.50 

2.2313E+06 

2.1044E-01 

4.9358E+00 

4 

2.00 

1 . 3534E+06 

2.3139E-01 

5.5195E+00 

5 

2.50 

8. 2085E+05 

1.8048E-01 

4.8660E+00 

6 

3.00 

4.9787E+05 

1.1483E-01 

4. 2230E+00 

7 

5.00 

6.7379E+04 

1.2439E-01 

2.2664E+00 

8 

8.00 

3.3546E+03 

7.2094E-03 

1.2042E+00 

9 

12.00 

6.1442E+01 

•0000E+00 

1.0533E+00 

10 

16.00 

1. 1254E+00 

•0000E+00 

1. OOOOE+OO 

11 

17.00 

4.1399E-01 

•0000E+00 

1.0266E+00 

12 

17.50 

2. 51 10E-01 

•0000E+00 

1. 1047E+00 

13 

18.00 

1.5230E-01 

.0000E+00 

1.4014E+00 

14 

18.50 

9.2374E-02 

.0000E+00 

1.8501E+00 

15 

19.00 

5.6028E-0 2 

•0000E+00 

1.8196E+00 

16 

19.50 

3.3983E-02 

•0000E+00 

1.2708E+00 

17 

20.00 

2.0612E-02 

•0000E+00 

6.8993E-01 

18 

20.50 

1.2502E-02 

•OOOOE+OO 

3» 1839E-01 

19 

21.00 

7.5826E-03 

•0000E+00 

1.3300E-01 

20 

21.50 

4.5991E-03 

•OOOOE+OO 

5.2364E-02 

AGE  TO  INDIUM  RESONANCE  U.46EV)  IS 
AGE  TO  ARBITRARY  THERMAL  (1.12EV)  IS 
THERMAL  DIFFUSION  LENGTH  SQUARED  IS 
TOTAL  MIGRATION  AREA  IS 

MOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS 
AVERAGE  THERMAL  NEUTRON  VELOCITY  IS 
ABSORPTION  PARAMETER  IS 
SCATTERING  PER  RESONANCE  ATOM  IS 
EFFECTIVE  RESONANCE  INTEGRAL  IS 


4.3711E+01  CM2 
4.3970E+01  CM2 
3.2857E+00  CM2 
4.7256E+01  CM2 
2.8887E+03  M/SEC 
3.8368E+03  M/SEC 
1.1892E+00 
6. F 078E+04  BARNS 
2.5021E+02  BARNS 


Figure  6a:  Output,  Chain  One 
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3 ■ 


100  1000 


NUCLIDE  IS  URANIUM235 


ITS  NUMBER  DENSITY  IS  1.164CE-04  PER  BARN-CM 


MACROSCOPIC  CROSS  SECTIONS,  IN  »ER  CM,  FOLLOW. 


GROUP 

NUSI6F 

SIGA 

SIGR 

SIGTR 

SOURCE 

1 

.000416 

.000159 

.000159 

.000157 

2.730272 

2 

.010673 

.005873 

.005873 

.001753 

,007209 

3 

.096155 

.045928 

.045928 

.001746 

.000000 

SCATTER 

FROM  TO 

INELASTIC 

N-2N 

ELASTIC 

CPO) 

ELASTIC 

(PI) 

1 

1 

.0000 

.0000 

.0000 

,0000 

1 

2 

.0000 

.0000 

.0000 

.0000 

1 

.0000 

.0000 

.0000 

.0000 

2 

2 

.0000 

.0000 

.0000 

.0000 

2 

3 

.0000 

.0000 

.0000 

.0000 

3 

3 

.0000 

.0000 

15.0000 

.0420 

NUCLIDE  IS  URANIUM  38 

ITS  NUMBER  DENSITY  IS  1.2900E-05  PER  3ARN-CM 
MACROSCOPIC  CROSS  SECTIONS,  IN  PER  CM,  FOt-LOW. 


GROUP 

NUSIGF 

SIGA 

SIGR 

SIGTR 

SOURCE 

1 

.000008 

.000004 

.000005 

.000073 

2.730272 

2 

.000000 

.000002 

.000002 

.000149 

.007209 

3 

.000000 

.000057 

.000057 

.000107 

.000000 

SCATTER 

FROM  TO 

INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 

(PI) 

1 

1 

.0000 

.0000 

.0001 

.0000 

1 

2 

.0000 

.0000 

.0000 

.0000 

1 

3 

.0000 

.0000 

.0000 

.0000 

2 

2 

.0000 

.0000 

.0002 

.0000 

2 

3 

.0000 

.0000 

.0000 

.0000 

3 

3 

.0000 

.0000 

8.3000 

.0232 

T 


LETHARGY 

5.000000 

17.000000 

99.999900 


TOTAL 

TRANSFER 

.0000 

.0000 

.0000 

.0000 

.0000 

15.0000 


LETHARGY 

5.000000 

17.000000 

99.999900 


TOTAL 

TRANSFER 

.0001 
.0000 
.OOOf. 
.0002 
.0000 
8. 3000 


NUCLIDE  IS  HYDROGEN 


ITS  NUMBER  DENSITY  IS  4.082CE-02  PER  BARN-CM 

MACROSCOPIC  CROSS  SECTIONS,  IN  PER  CM,  FOLLOW. 

GROUP  NUSIGF  SIGA  SIGR  SIGTR 

SO,,r  ” 

LETHARGY 

1 

.000000  .000000 

•  054383 

.050234 

2  730272 

5.000000 

2 

.000000  .001101 

.208083 

.273940 

.007209 

17.000000 

3 

.000000  .007714 

.007714 

.525538 

.000000 

99.999900 

SCATTER 
FROM  TO 

INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 

(PI) 

TOTAL 

TRANSFER 

1  1 

.0000 

.0000 

.  1552 

.1152 

.1552 

1  2 

.0000 

.0000 

.0544 

.0231 

.0544 

1  3 

.0000 

.0000 

.0000 

.0000 

.5000 

2  2 

.0000 

.0000 

.6116 

.4365 

.6116 

2  3 

.0000 

.0000 

.2070 

.1081 

.2070 

3  3 

.0000 

.0000 

38.0000 

25.1332 

38.0000 

NUCLIDE  IS  OXYGEN 

ITS  NUMBER  DENSITY  IS  2.4100E-02  PER  BARN-CM 

MACROSCOPIC  CROSS  SECTIONS,  IN  PER  CM,  FOLLOW. 

GROUP  NUSIGF  SIGA  SIGR  SIGTR 

SOURCE 

LETHARGY 

1 

.000000  .000206 

.001271 

.054872 

2.730272 

5.000000 

2 

.000000  .000000 

.003552 

.084439 

.007209 

17.000000 

3 

.000000  .000003 

.000003 

.096959 

.000000 

99.999900 

SCATTER 
FROM  TO 

INELASTIC 

N-2N 

ELASTIC 

(PO) 

ELASTIC 

(PI) 

TOTAL 

TRANSFER 

1  I 

.0000 

.0000 

.0766 

.0080 

.0766 

1  2 

.  0000 

.0000 

.0011 

-.0003 

.0011 

1  3 

.0000 

.0000 

.0000 

.0000 

.0000 

2  2 

.0000 

.0000 

.0846 

e  OL  47 

.0846 

2  3 

.0000 

.0000 

.0036 

-.0011 

.0036 

3  3 

.0000 

.0000 

4.2000 

.1751 

4.2000 

Figure  6d:  Output,  Chain  Two,  Continued 
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NUCLIDE  IS  ALUMINUM 


ITS  NUMBER  DENSITY  IS  2.3330E-02  PER  BARN-CM 


MACROSCOPIC  CROSS  SECTIONS.  IN  PER  CM,  FOLLOW. 


GROUP 

NUSIGF  SIGA 

SIGR 

SIGTR 

SOURCE 

1 

.000000  .000073 

.001230 

.057909 

2.730272 

2 

.000000  .000527 

.001230 

.032878 

.007209 

3 

.000000  .003073 

.003073 

.031525 

.000000 

SCATTER 

ELASTIC 

ELASTIC 

FROM  7 

0 

INELASTIC 

N-2N 

(PO) 

(PI) 

1 

1 

.0044 

.0000 

.0768 

.0184 

1 

2 

.0000 

.0000 

.0012 

-.0004 

1 

3 

.0000 

.0000 

.0000 

.0000 

2 

2 

.0000 

.0000 

.0325 

*0010 

2 

3 

.0000 

.0000 

.0007 

-.0002 

3 

3 

"OCCO 

.0000 

1.4000 

.0518 

MAXWELL 

-BOLTZMAN  FACTOR  =  1.123 

,  #4HHW-WIGNER  =  1.7462097E+00 

Figure  6e:  Output,  Chain  Two,  Continued 


i 

t 

i 


i 


I 


lethargy 

5.000000 

17.000000 

99.999900 


TOTAL 

TRANSFER 

.0812 
.0012 
•  0000 
.0325 
.0007 
1.4000 
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APPENDIX  A 

JOB  OLD  BARNYARD  CHAIN  1*  MAY  1966 

THIS  CARD  ORIGINS  THE  PROGRAM  AT  ADDRESS  10970  IN  MEMORY.  IN  ORDER 
TO  DO  THIS  YOU  MUST  USE  A  SPECIAL  SET  OF  SHORT  SUBROUTINES. 

ORIGIN  10970 

STORING  IN  COMMON  THOSE  DATA  USED  IN  CHAIN  2  OR  IN  SUBPROGRAMS 
COMMON  LGP( 13) .KE.FLX! 11 ) ,U! 11 ) * SOURCE! 11 ) .DENaR .NRNUC .TEMPTR  »XBAR 
DIMENSION  SA (11) .STR! 11 ) ♦SR! 11) »XSAM(168) »SSN(12»13) ,S2N(12,13) »X( 

1 1 1 )  ♦  T  ( 1 2  ) 

PEADIN6  IN  THE  FIRST  TWO  DATA  CARDS 
I EAD  1060. LGP 
I  FORMAT! 1 3A5 ) 

READ .NEW .NRNUC .DENBR . TEMP TR .RES 
PUNCH1060 . LGP 
PUNCH1065 
FORMAT l//) 

TYPE1060.LGP 

READING  IN  THE  DECK  OF  SOURCES 
READ  1070. U 
K  =  0 

DO  1046  J=1.10 
IF'K'1045. 1043, 1045 
.  READ1060.LGP 
READ1070, INDEX, SOURCE 
i  FORMAT! 12N) 

I F  C 1 NDEX-NEW ) 1 046 , 1 044 , 1 046 
K=1 

GO  TO  1046 
READ  1060, PU 
READ.PU 
CONTINUE 

I F(NEW-IC) 1050,1049.1047 

SETTING  THE  UNIT  SOURCE  IN  ITS  GROUP  (IF  APPLICABLE) 

SOURCE (NEW-10)=1. 

LGP! 10)=NEW-10 
PUNCH  1071 .LGP 
FORMAT (9A5.I3.3A5) 

GO  TO  1048 

READING  IN  THE  USERS  OWN  SOURCE  (IF  APPLICABLE) 

READ107Q. SOURCE 
PUNCH  1060 »LGP 
TYPE  600 

FORMAT (/16HBEG IN  EPITHERMAL) 

PUNCH  1058.TEMPTR 

FORMAT! /22HSYSTEM  TEMPERATURE  IS  .F6.1.7H  KELVIN) 

PUNCH  1059 

i  ORMAT! /18HAT0M  FRACTIONS  ARE/) 

CLEARING  THE  MATRICES  OF  ANY  PREVIOUS  PROBLEM 
DO  1000  J=l,ll 
STR ( J ) =0. 

SR! J)=0. 


'[I  WJ! 


a 


\  '  l 


■\* 


1000  SR ( J  j - • 0 

DO  1002  J= 1  *  12 
DO  1002  K= 1 > 1 3 
SSN( J,K)=0. 
1002  S2N (  J  *K ) =0, 


I  READING  IN  THE  RESONANCE  DATA  CARD  (IF  APPLICABLE) 

IF!RES)995, 994,995 
995  READ*JJJ*PU»RIEFF 
PU=PU/DEN3R 


SETTING  UP  THE  DO  LOOP  FOR  EACH  OF  THE  NUCLIDES 
PROB=l • 

KKK=0 

PSIGF=0. 

PSIGM=0. 

PSIG0=0. 

SIGTT=.0 
SIGDEL=«  0 
SIGSTH=0. 

DO  1020  1  =  1  »NRN'JC 


I EADING  IN  THE  NUCLIDE  DATA  CARDS 
READ  1061*  SAM,  SAMM,  ATOMS,  P,  MOD,  DENS 
1061  FORMAT(3A5,9N) 

P=P/DENBR 

READ  1061, SAM, SAMM, ATOMS, LS, NA, LQ, PSI ,SIGT,SIGA, DELTA, SIGS 
PUNCH 10 66, SAM, SAMM, ATOMS, P 
1066  FORMAT (3A5,F15»8) 


CALCULATING  TOTAL  THERMAL  SCATTERING  SIGMA 
SIGSTH=SIGSTH+P*SIGS 


:  CALCULATING  THERMAL  UTILIZATION  AND  ABSORPTION  PARAMETER 

GO  TO  (993,992,996) , MOD 

992  PSIGF=PSIGF+P*SIGA 
GO  TO  991 

996  PSIGO=PSlGO+P*SIGA 
GO  TO  991 

993  PSIGM=PSIGM+P*SIGA 
SIGDEL=SIGDEL+SIGT*P*DELTA 
SIGTT=SIGTT+P*SIGT 


C  READING  IN  THE  WHOLE  CROSS  SECTION  LIBRARY  FOR  THE  NUCLIDE 

991  READ  1067, (XSAM( J) , J=1,LS) 

1067  FORMAT (6N) 

KZ  =  0 
KKZ  =  0 


CHECKING  FOR  ABSOR°TION 
IF(NA)1006. 1006, 1004 

1004  DO1005  J=(,ll 
KZ  =KZ  + 1 

STR( J)=STR( J)4XSAMCKZ)*P 

1005  SR  (  J )  -  S R  (  J  )  +  X  $ AM ( KZ ) * P 


3U 


1006  KZ=KZ+1 

I N=XSAM ( KZ ) 

KZ- KZ+1 
I  IN  =  XSAM(KZ ) 

DO  1017  J=1 » IN 
KK=J+I IN 

I F  (  KK-  x i.  1 1008  *  1008  * 1007 

1007  KK  = 1 2 

1008  DO  1C  1 7  K  =  J*<K 

GO  TO  « 1009,1010*1009,1011 ) *L0 

1009  KZ=KZ+ 1 

C  INELASTIC  CROSS  SECTION*  GROUP  TO  GROUP 

SSN ( J ,K+1 ) =SSN (J»K+1)+XSAM(KZ)#P 
GO  TO  1012 

1010  KZ=KZ+1 

C  N-2N  CROSS  SECTIONS*  GROUP  TO  GROUP 

S2N( J*K+1 )=S2N(J»K+1)+XSAM(KZ)*P 
GO  TO  1012 

1011  IF(KKZ)2001*2000»2001 

2000  KZ=KZ+1 
KKZ=  1 
GOT 02002 

2001  KZ=KZ+2 

C  PUTTING  ELASTIC  (PO)  AND  ELASTIC  (PI)  IN  THE  BOTTOM  OF  SSN  AND  S2N 

2002  LS= 13— J 
LE= 13-K 

SSN(LS»LE)=SSN(LS,LE)+XSAM(KZ)*P 

S2N(LS»LE)=S2N(LS*LE)+XSAM(KZ+l)*P/3. 

1012  STR( J)=STR( J)+XSAM(KZ)*P 
IF(K-J)1016, 1013, 1016 

1013  r-0  TO  (1017, 1014*1017*1015), LQ 

1014  .*.R(  J)=SR(  J)-XSAM(KZ)*P 
GO  TO  1017 

1015  STR( J)=STR( J)-XSAM( KZ+l )*P/3. 

GO  TO  1017 

1016  SR( J)=SR( J)+XSAM(KZ )»P 

1017  CONTINUE 

GO  70(1018*1018,1019,700) *LQ 

1018  LQ=4 

GO  TO  1006 

1019  LQ=2 

GO  TO  1006 

C  READING  IN  THE  GRID  BACKGROUND  FOR  THE  SPECTRUM  PRINTOUT 

700  READ  2222* (XSAM(J) »J=41, 120) 

2222  FORMAT ( 80A1 ) 

C  CHECKING  FOR  RESONANCE  CALCULATION 

IF(KKK)705, 701*705 

701  IF (RES) 704*1020,704 
704  KKK= 1 
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C  CALCULATING  PSI*SIGMA-S  FOR  THE  MIXTURE 

TS I G8=0 • 

TSIG9=0. 

P’S  I  G8  =  0  • 

705  SIG8=SSN(5,5)+SSN(5»4)+SSN(5,3)+SSN( 5.2)+SSN(5»l )-TSIG8 

TSIG8=SSN(5»5)+SSN( 5,4)+SSN(5  ,3)+SSN( 5,2>+SSN(5»l ) 

SIG9  =  SSN(4,4)+SSN( 4,3)+SSN<4,2  >+SSN(4, 1 )~TSIG9 
TSIG9=SSN(4,4)+SSN( 4,3)+SSN(4 ,2 )+SSNC4, 1 ) 
PSIG8=PSIG8+PS1*P*(SIG8+SIG9)*.5 
1020  CONTINUE 

C  ALL  CROSS  SECTIONS  HAVE  BEEN  CALCULATED  AT  THIS  POINT 

IF(KKK) 720,730,720 

720  TYPE  601 

601  FORMAT ( /15H8EGIN  RESONANCE) 

C  CALCULATING  SCATTERING  PER  RESONANCE  ATOM  (IF  APPLICABLE) 

SCAT=.5*(TSIG8+TSiG9) /PU 

C  SELECT ’ NG  THE  EMPIRICAL  EQUATION  FOR  RESONANCE  INTEGRAL 

IF  <RIEFF)728*721,728 

721  GO  TO(722*?25) »JJJ 

722  IFCSCAT-4000. >723,723,724 

723  RIEFF=2.69*SCAT**.471 
GO  TO  728 

724  RIEFF  =  EXPF(5»64-163«/(SCAT*>f»65)  , 

GO  TO  728 

725  IF(SCAT-4500. 1726.726,727 

726  RIEFF=8.33«SCAT**.253 
GO  TO  728 

727  R I EFF- 70 • 

C  CALCULATING  RESONANCE  ESCAPE  PROBABILITY 

728  PR08=EXPF(PU*RirFF/PSIG8 ) 

C  ADJUSTING  GROUP  8  REMOVAL  CROSS  SECTION  FOR  RESONANCE  ESCAPE 

SR  1 8 ) =SR ( 8 ) *  PROB 
STR(8)=SR(S;+S5N(5,5)-S2N<5,5) 

C  HANSENS  FINITE  DIFFERENCE  EQUATIONS  FOR  CPI  THERMAL  FLUX  AND  AGE 

73C  SAMM=1 ./  (3.*Dl'NBR*DENPR  ) 

PUNCH  1051 

1051  FORMAT  (//57X.14HRELATIVE  GROUP/ 5HGROUP ,  7X  » 8HLETHARGY  ,  7X  ,  lOHENER'oY  , 
1  EV,8X,6HS0URCE  ,6X,  14HFLUX,  /UNIT  U'/> 

RHO=.C 

DO  1029  K  = 1 ,11 
LE= 1 3-K 

DENS=SAMM/(STR(K)*SR(K) ) 

AT0MS=C. 

TOTWT  =0. 

P  =  0. 

T ( K ) =0» 

FLX ( K ) =C • 

SMK)=0. 
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IF ( K  —  1 ) 1027,1027,1025 
10P5  DO  1026  J- 1  * K- 1 
L  S= 1 3- J 

SAM=TLX  (  J )  *  {  5S.N  (  L S , L E  )  +SSN  <  J ,.<+ 1 )  +2  « *S2N  (  J,K U  ) )  /SR  (  J  ) 

AT0M3=AT0MS+SAM 

TOTWT=TOTWT+SAM*T ( J ) 

1026  P=P+S2N(LS»LE)*SA( J) 

1027  FLX ( K ) = A TOMS i SOURCE ( K ) 

3A t  K )  =  FLX ( K ) *DENS*P/STR ( K ) 

T ( K ) = ( SA ( K ) +TOTWT ) /FLX ( K ) 

X ( K ) =U( K ) -RHO 
RHO=U ( K ) 

1029  CONTINUE 

F=SR( 10 ) *4 • /FLX ( 10) 

C  CALCULATING  TOTAL  FLUX  PER  GROUP 

DO  1023  K=1 ,10 

1023  FLX(K)=FlX(K)/SR(K) 

C  CALCULATING  AND  NORMALIZING  FLUX  PER  UNIT  LETHARGY 

DO  1022  K=1 » 10 
XSAM ( 40  ~K ) =F*FLX ( K ) /X(K) 

KK=40-K 

C  OUTPUTTING  EPITHERMAL  RESULTS 

1022  PUNCH  1052tX»U(K)»1.0E+07*EXP<-U(K)>.SOURCE(K>»XSAM(KK) 

C  NOTICE  THAT  TOTAL  FLUX  PER  GROUP  HAS  BEEN  LEFT  IN  STORAGE  IN  FLX(ll) 

C  CALCULATING  AGE  TO  INDIUM  RESONANCE 

SAMM=(1  ( 10  ) -T  (  9 )  )*0«  935-»-T  (  9  J 

C  (HECKING  FOR  MODERATOR  IN  THE  SYSTEM 

I F(PSIGM) 1024, 1021, 1024 
1021  PUNCH  850 , SAMM 
TYPE  1028 

1028  FORMAT (36HTHERE  IS  NO  MODERATOR  IN  THE  SYSTEM. /36HN0  THERMAL  CALCU 
1LATI0N  WILL  BE  DONE.) 

GO  TO  860 

C  STARTING  THE  THERMAL  CALCULATION 

1024  TYPE  602 

602  FORMAT ( /13H8E6IN  THERMAL) 

C  THE  ABSORPTION  PARAMETER 

DELTA=SIGDEL*(PSIGF+PSIGM+PSIG0)/(S1GTT*(P5IGM+PSIG0) ) 

AREAX=.0 
PREFX-.O 
PREVF=0 • 

11=0 

AN2=1. 

P  =  0. 

C  ARRANGING  FOR  INCREMENTS  OF  0.25  LETHARGY  UNITS 
ETODU=EXP( .125 )-l. 
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V  HERE  EQUAL F  Xo  clNCF  X  WAS  SUB  SC  R  i  P  T  i  D  ABOVE 

Y  =  E  XP ( -  1 ] • S } *  SOR  T  ;  1 .OF+i?/ ( 8. 6 16* TEMPI R )  ) 

1=0 

IOING  THE  SERIES  SOLUTION  TO  WlLKiNS  EQUATION 

AN1=AN2*DELT A*. 333333 

DSUM  =  AN  1 

SUM=AN2+AN1*Y 

XTON=Y 

AN3=0. 

N  =  4 
N  =  N-i  1 
XN-N-3 

IF  ( P) 805 >804,805 

THE  RECURSION  RELATION  FOR  X=0 

AN=  I  2.*(XN-2.  »  *AN?+DELTA*AN1  ) /( XN*  (  XN  +  2.  )  ) 

GO  TO  606 

THE  GENERAL  RECURSION  RELATION 

AN* I  AMI  * ( XN-1 • )*<XN+l.-2.*P*P)-AN2*(4.*P*(XN-2» ) +0ELTA ) -AN3*2 •  * ( X 
IN-3. ) ) /vP*XN*< l.-XN) ) 

AN3= AN2 

DTERM=XN*AN*XTON 

XION=XTON*Y 

TERM=AN*XTON 

SUM=SUM+TERM 

DSUM=DTERM+DSUM 

AN2=AN1 

AN1=AN 

CHECKING  FOR  CONVERGENCE  OF  THE  SERIES 
IF ( ABS ( TERM/SUM) -1.0E-06 5 807* 807, 803 
IF(ABS(DTERM/DSUM I- 1.0E-06) 808,808,803 
P  =  Y+P 

FINDING  N ( X )  FOR  THE  POINT  X=P 
F=EXP ( -P*P ) *SUM*P**2 

LOCATING  THE  PEAK  OF  THE  NIX)  CURVE 
I F ( I  I >816,814,816 
IF(F-PREVF)315,815,816 

(ALCULATING  THE  MOST  PROBABLE  VELOCITY 
VPROB  =  SQR  T ( 1.6477E+04*TEMPTR)*P 
11  =  1 
PREVF  =  F 

CALCULATING  FLU.X  PER  UNIT  U  WHICH  IS  X*X*N(X) 

F=F*P*P 
1  =  1  +  1 

STORING  THE  VALUES  OF  THERMAL  FLUX 
XSAMI I ) =F 
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C  CALCULATING  THE  area  UNDER  THE  XN  (  X  )  CURVE  BY  TRAPEZOIDAL  RULE 

AREAX=AREAX+(F+PREFX) *Y*.5 
PREFX=F 

IFIP-6. 6)809, 810*810 

809  AN2=SUM 
AM 1 =DSUM 
Y=P+P*ET0DU 
Y  =  Y-P 

GO  TO  802 

C  CALCULATING  THE  AREA  UNDER  THE  FLUX  CURVE  ANALYTICALLY 

810  C=2.*F»( l./(P*P)-l. )+p*p*EXP!-P*P)*DSUM 
AREA=(P*C+2.*(P*P-1. )*F/P) /DELTA 

C  FINDING  THE  AVERAGE  NORMALIZED  VELOCITY 

;.BAR=AREAX/AR£A 

C  NORMALIZING  THE  STORED  THERMAL  FLUXES 

F=1./XSAM(29) 

DO  312  J=l,29 
812  XSAM( JI=XSAM( J)*F 
P=16.5 


C  OUTPUTTING  THE  THERMAL  FLUXES 

DO  813  K=1 ♦ 10 
P  =  P+ • 5 
J=2 7-K-K 

813  PUNCH  1052,K+10*P» 1.0E+07*EXP(-P) .0® »XSAM! J) 

1052  FORMAT! 1 3 » FI  5. 2 , 3X , 1P3E16. 4 ) 

C  OUTPUTTING  THE  THERMAL  CONSTANTS 

PUNCH  850  *SAMM 
PUNCH  85 1 »  T ( 10 ) 

850  FORMAT (/ /5X»35HAGE  TO  INDIUM  RESONANCE  (1.46EV)  IS » 1PE20.4 ,4H  CM2) 

851  FORMAT! /5X*36HAGE  TO  ARBITRARY  THERMAL  !1.12EV)  IS*  1PF19,4,&h  CM2 
1) 


C  CALCULATING  THE  DIFFUSION  LENGTH  SQUARED  AND  MIGRATION  AREA 

S I GATH=PS I GM+PS I GF+PS I  GO 

ELSQRD=. 33333*SIGSTH*XBAR/! SI  GATH* !DENBR*(SIGSTH+SIGATH/XBAR) }**2) 

PUNCH  852 tELSORD 

852  FORMAT! /5X , 3 5HTHERMAL  DIFFUSION  LENGTH  SQUARED  I S» 1PE20.4 ,4H  CM2) 
EMSQRD=T ( 10 ) +ELSQRD 

PUNCH  853 » EMSQRD 

853  FORMAT! /5X,23HTOTAL  MIGRATION  AREA  IS* 1PE32.4,4H  CM2) 

PUNCH  854 ♦ VPROB 

854  FORMAT! /5X»41HMOST  PROBABLE  THERMAL  NEUTRON  VELOCITY  IS* 1PE14.4.6H 
1  M/SEC) 

VBAR  =SQRT ! 1 .6477E+04*TEMPTR ) *X3AR 
PUNCH  855.V8AR 

855  FORMAT! / 5X , 3 5HAVERAGE  THERMAL  NEUTRON  VELOCITY  IS* 1PE2C*4*6H  M/SEC 
1) 

PUNCH  856>DELTA 

FORMAT! /5X,23HABSORPTION  PARAMETER  IS»1PE32.4) 
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C  OUTPUTTING  THE  RESULTS  Or  THE  RESONANCE  CALCULATION 

IF'PP S)859, 860,859 
HS9  PUNlh  8  5  7  »  SC  A  T 

857  F  ORMAT I /5X  *  32HSCAT TER  I  MG  PER  RESONANCE  ATOM  I S » 1 PE2 3  *4  » 6H  BARNS) 
PUNCH  858.RIEFF 

858  FORMAT < /5X»31HEFFECT I VF  RESONANCE  INTEGRAL  I S» 1PE24.4  »6H  BARNS) 

860  PUNCH  861 

861  FORMAT ! 78X  » 1H- ) 

C  WARNING  THE  USER  TO  SET  SWITCH  FOR  SPECTRUM  OUTPUT 

TYPE  362 

862  FORMAT ( /52HPUT  SW’TCH  1  ON  FOR  SPECTkUM  PRINTOUT.  PRESS  START.) 
PAUSE 

I F ( SENSE  SWITCH  1)3000.2909 

C  THE  SPECTRUM  PUNCHOUT  SUBPROGRAM 

2999  TYPE  863 

863  FORMAT! /28HNO  SPECTRUM  WiLL  5E  PUNCHED.) 

GO  TO  1091 

3000  TYPE  603 

603  FORMAT! /14HSIGIN  SPECTRUM) 

C  STORING  THE  SYMBOLS  FOR  THE  GRAPH 

X(1)=1HX 
X!2)=1HL 
X ( 3 ) =  1HH 
PUNCH  300 

300  ;  ORMAT (69HPL0T  OF  SPECTRUM,  IN  FLUX  PER  UNIT  LETHARGY  VERSUS  LETHA 

1RGY  *  FOLLOWS.//) 

PUNCH  301 

C  REDUCING  FLUXES  TO  LOG  SCALE  80  UNITS  WIDE 
DO  jOOI  J=)  ,39 

3C01  XSAM(J)  =  f ! LOG <  XSAM ! J ) ! /2 . 30259 )  + 1 . ) *20 .  +  . 5 
KK  =  9 
1=0 

C  ONE  PASS  THROUGH  THE  LOOP  FOR  EACH  LINE  OF  OUTPUT 

DO  3005  J=1 ,93 
IF(KK)30U, 3014, 3011 

C  CHECKING  FOR  EPITHERMAL  OR  THERMAL  FLUXES 

3011  IF! J-29J3009, 3009, 3012 

3012  PU=J 

P  =23.2-PU/4. 

C  ADJUSTING  FOR  DIFFERENT  WIDTH  GROUPS 

I F { P-U ( KK ) ) 301 3, 30 13, 301 4 
3C1 3  KK=KK-1 
301 A  K  =  XSAM  t  39-KK ) 

GO  TO  3008 
3009  K=XSAM(J) 

3008  1=1+1 


NUMBERS  FOR  THE  LETHARGY  SCALE 
1 E-23-J/4 

CHECKING  FOR  VALUES  TOO  LARGE  FOR  SCALE 
IF(K-79)3015, 3016, 3006 
3006  N  =  3 
<  =  79 

GO  TO  3004 

CHECKING  FOR  VALUES  TOO  SMALL  FOR  SCALE 

3015  IF(K-3) 3019,3017,3016 

3017  GO  TO  (3018,3016,3016,3016,30181,1 

3018  N=2 
K  =  4 

GO  TO  3004 

3019  IF(K-2) 3018,3017,3017 

3016  N= 1 

3004  K2=K+39 
KK2  =  K  +4 1 

GO  TO (3002, 3003, 3003,3003, 3002) ,1 
OUTPUTTING  THE  SPECTRUM  GRAPH 

3  0C2  PUNCH  302, LE, (XSAM(LS) ,LS=4  3,KZ) *  X ( N  » , (XSAM(LS) ,LS=KKZ»120) 

302  FORMAT! I2,78A1 ) 

1  =  1 

GO  TO  3003 

3003  PUNCH  2222, (XSAM(LS) ,LS=41,KZ) ,X(N) * ( XSAM ( LS ) ,LS=KKZ » 120 ) 

3005  CONTINUE 
PUNCH  301 

301  FORMA"!  /2H.1,17X,1H1,18X»2H10,18X,3H100,15X»4H1000/) 

1091  TYPE  1092 

1092  FORMAT ( /14HEND  OF  CHAIN  1//39HVC  COuLAPSE  CROSS  SECTIONS  LOAD  CHAI 
IN  2/1 

THIS  PUNCH  MAKES  THE  407  LISTER  ADVANCE  A  SHEET  OF  PAPER 
PUNCH  861 

THIS  PUNCH  AVOIDS  RUNNING  OUT  THE  LAST  PUNCHED  OUTPUT  CARD 
PUNCH  1065 


THIS  IS  DAR T  II  OF  OLD  BARNYARD.  CROSS  SECTION  COLLAPSING 


S  JOB  OLD  BARNYARD  CHAIN  2 

T.  ORIGIN  10970 

C  NOTICE  COMMON  STATEMENT  IS  SAME  AS  IN  CHAIN  1 

COMMON  LGP( 13) *KE * FLX (11) ,UC 11 ) .SOURCE (11) .DENBR ,NRNUC ♦ TEMPTR ,XBAR 
DIMENSION  SA{ 11 ) ,STR( 1 1 >  *  SR  Ill? .XSAM( 168) ,SSN( 12.13) ,SZNC 12.13), 
1SNC 12,13) 

TYPE  604 

604  FORMAT ( /16H8EGIN  COLLAPSING) 

C  READING  IN  NUMBFR  OF  BROAD  GROUPS  AND  THEIR  BOUNDARY  FINE  GROUPS 

READ  1070.NBG, (LGP( J) ,J=2,NBG+1) 

1070  FORMAT ( 12N ) 

MXL=NBG+1 
LGP ( 1 ) =0 
LGP ( N3G+2 ) = 12 

C  BETTER  LOOK  NOW  AT  FUNCTION  SUBROUTINE  KS  AT  END  OF  PROGRAM 

-  NORMALIZING  THE  FINE  GROUP  FLUXES  IN  EACH  BROAD  GROUP 

:  AND  ADDING  UP  THE  SOURCE  IN  EACH  BROAD  GROUP 

DO  2  J= 1 »NBG 
SAC J)=0. 

DO  2  K=KS( J) ,KE 

2  rAiJ'=SA(J)+FLX(K) 

I  0  4  J= 1  ,N3G 

SAM= • 0 

DO  3  K=KS ( J ) ,KE 
SAM=SAM+SOURCECK) 

3  CLXCK)=FLX(K)/SA(J] 

4  SOURCE ( J)=SAM 

:  ONE  TRIP  THORUGH  THE  DO  LOOP  FOR  EACH  NUCLIDE 

DO  55  N= 1 .NRNUC 

l  CLEARING  THE  MATRICES  OF  ANY  PREVIOUS  PROBLEM 

DO  555  J=l,ll 
SAC J)=0. 

SRC J)=C. 

555  STR ( J ) =0. 

DO  6  J=i  ,  12 
DC  6  K  =  1 , 13 
SSNC J.K) =0r 
6  S2N ( J  »K ) =0. 


READING  IN  THE  NUCLIDE  DATA  CARDS 
READ  106 1, SAM, SAMM, ATOMS ,P, MOD .DENS 
1061  FORMAT ( 3A5 ,9N ) 

P=P/DENBR 

READ  1061,SAM,SAMM,AlOHS,LS,NA,LG,PSI ,SIGT,SIGA,DELTA,SIGS»T£RMU 
DUNCH  200, SAMM, ATOMS, P*DENBR 

200  FORMAT { //l 1HNUCL IDE  IS  ,2A5,///25H  ITS  NUMBER  DENSITY  IS,1PE11. 
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l'-.^H  PER  BARN-CM//} 

C  CHECKING  FOR  MICRO  OR  MACROSCOPIC  CROSS  SECTIONS 

I F ( DENS-i • 161.62,61 
t.1  DENS=P*DENBR 
PUNCH  20? 

202  FORMAT { 4S-HMACR0SC0P I C  CROSS  SECTIONS.  IN  PER  C.M.  FOLLOW..) 

GO  TO  49 
62  PUNCH  201 

201  FORMAT ( 45HMI CROSCOPI C  CROSS  SECTIONS.  IN  BARNS.  FOLLOW.) 

49  PUNCH  220 

220  FORMAT ( //79HGROUP  ABSORPTION  REMOVAL  TRANSPORT  DlFF  CO 

1EF  SOURCE  LETHARGY/) 

C  READING  IN  THE  WHOLE  CROSS  SECTION  LIBRARY  FOR  THE  NUCLIDE 

READ  1067. IXSAMI J; .  '  '  'A) 

1067  FORMATION) 

KZ  =  0 

C  CHECKING  FOR  ABSORPTION 

IFCNA) 10.10.7 
7  DO  9  J=1,N3G 

DO  3  K=KS; J) .KE 
K?  =KZt-l 

SA  (  J  )  =SA  (  J  )  +XSAM  (  KZ  )  *FLX  (  K  )  *Dt‘NS 
S  STR ( K  5  =XSAM '  KZ  ) 

9  SR ( J ) =SA ( J ) 

10  •  DO  12  J=  1 . 12 

DO  1?  K= 1 » 13 

12  SN{ J.K) =0» 

KZ=KZ+1 

IN=XSAM(KZ) 

KZ=KZ+1 
IIN=XSAM(KZ5 
DO  15  J= 1 . 1 N 
KK=J+I IN 

1 F ( KK-1 2 )14, 14.13 

13  KK= 12 

14  DO  15  K  =  J » KK 

r-0  TO  (  141,141,141,371)  ,LO 

C  USING  SN  AS- A  TEMPORARY  STORAGE  MATRIX  FOR  INELASTIC  OR  N-2N 

141  KZ-KZ+1 

STRI J)=STR{ J)+XSAM(KZ! 

SN( J.K+1 ) =XSAM ( KZ J *FI  X( J)*DENS 
GO  TO  15 

C  USING  SN  AS  A  TEMPORARY  STORAGE  MATRIX  FOR  ELASTIC  <P0)  AND  (PI) 

371  NA=13-J 

KE=13-K 
KZ=KZ+1 

XSAMIKZ  +  1)=XSAM(KZJ-1)*«33333 
STR  f  J )  =  STR  (  J  )+-XSAM(  KZ )  -XSAMf  KZ+1 ) 

SN ( J.K+1 I -XSAM ( KZ ) *FLX ( J ) *D£NS 
SN(NA,KE)-XSAM(KZ+1)*-FLX(J)*DENS 
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K2=K?.  +  1 
15  CONTINUE 

DO  19  J=?.,NB6 
DO  19  <= J»NBG+1 
LS-LCP ( X )  +1 
L£-«.GD(  K+ 1 » 

DO  19  i=KS(  I)  »£& 

I  F { LS-I 5 16,17,17 
LS=  I 

DO  19  L  =  LS »LE 
GO  TO  ( 11,22*11,34) »L0 

T NELAST I C  CROSS  SECTION ,  GROUP  TO  GROUP 
i SN( J,X+1 ) =SSN ( J  »K+1 ) +SN ( I » L+ 1 ) 

GO  TO  99 

ELASTIC  CROSS  SECTION*  GROUP  TO  GROUP  ( s>0 )  AND  (Pi) 

NA= 13— J 
KK=13— K 
IN=13-I 
1 1 N= 13-L 

SSN( NA*KK ) =SSN (NA  »KK ) +SN ! I  *  L+l ) 

S2N(NA,KK)=S2N(NA,KK)+SN{ IN, I  IN) 

99  IFCJ-K) 10,19,19 

I  N-2N  CROSS  SECTION,  GROUP  TO  GROUP 

22  32N( J ,K+1 ')  -S2N  ( J  »K+1 5+SN  ( I ,  L+l  5 

I F ( J— K) 18,30*18 
30  SR  ( J )  =SR  ( J )  ~SN  C I  »  L-*  1 ) 

GO  TO  19 

18  SR(J)=SR(JU-SNU,U1) 

19  CONTINUE 

GO  TO  (91,91,92,421) , LG 
9]  uQ=4 

GO  TO  10 
92  LG=2 

GO  TO  10 

Z  CALCULATING  BROAD  GROUP  SIGMA  TRANSPORT 

421  DO  46  J=1,NBG 

S4M=0. 

DO  45  K=KS( J) »KE 

Z  NOTICE  RECIPROCAL  RELATIONSHIP 

4;  SAM=SAM+FLX(K)/STR(K) 

46  3TR( J)=DENS/3AM 

3  OUTPUTTING  SINGLE  GROUP  VALUES 

DC  50  J=1*N8G 
K  =  U+1 

50  PUNCH  102  *  J  ,SA  ( J  )  ,SR(J)  ,STR  (  J )  ,  ,  33333/STR  ( J )  , SOURCE!  J)  »U(LGP(K)  *, 

I  CALCULATING  AND  OUTPUTTING  THERMAL  GROUP  VALUES 

SIGA=SI6A/X8AR 
SIGTR=SIGS*TERMU+SIGA 

4li 


PUNCH  102*MXL,SIGA,SIGA,S1GTR,.33333/SIGTR>0*»99.9999 

102  FORMAT ( I3»3X»6F12,4) 

C  READING  AND  IGNORING  THE  GRID  BACKGROUND  CARD  IN  EACH  library  DECK 

READ  1061, SAM 
PUNCH  205 

2C5  FORMAT ( / /7H3CAT  f ER »  37X  » 7HEI.AST i C» 7X*7HELASTI C , 8X  »5HTOTAL/8HFROM  T 
1C»7X,<JH  INELASTIC*  8X  »4HN-2N » 10X  *4H<  PO  i  »10X,4H(P1,  ,  7X»8HTRANSFER/  ) 

C  OUTPUTTING  GROUP  TO  GROUP  VALUES 

DO  54  J=1,NBG 
DO  5'*  K=J»MXL 
LS-13-J 
LE=13-K 

SAM  =SSN( JfK+1 )+2.*S2N( J ,K+1 ) +SSN ( LS » LE ) 

54  PUNCH103*J»K*  SSN ( J » K  +  3 ) ,S2N! J,K+1) » SSN ! LS , LE ) » S2N ( LS , LE ) »SAM 

C  OUTPUTTING  THERMAL  GROUP  TO  GROUP  VALUES 

55  PUNCH  103»MXL»MXL»0. *0. »SIGS» ( 1 » -TERMU ) *SI GS » SiGS 

103  FORMAT! 13,15, 5F14. 4) 

1093  TYPE  1063 

1C63  FORMAT! /15HEND  OF  PROGRAM,/) 

PUNCH  1065 

1065  FORMAT! //IX) 

STOP 
i  ND 

C  A  FUNCTION  SUBROUTINE  TO  SAVE  A  LOT  OF  WRITING  ON  DO  LOOPS 

FUNCTION  KS(J) 

COMMON  LGP ! 1 3 ) ,KE 
KS=LGP ( J ) +1 
KE=LGP( J+l ) 

RETURN 

END 
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